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(57) ABSTRACT 

A system and method for analyzing a surface. The system 
includes a computer including a CPU and a memory 
medium operable to store programs executable by the CPU 
to perform the method. The method may include: 1) receiv- 
ing data describing an n-dimensional surface defined in a 
bounded n-dimensional space, where the surface is embed- 
ded in an m-dimensional real space via embedding function 
x( ), and where m>n; 2) determining a diffeomorphism f of 
the n-dimensional space; 3) computing the inverse transform 
r 3 of the diffeomorphism f; 4) selecting points, e.g., a Low 
Discrepancy Sequence, in the n-dimensional space; 5) map- 
ping the points onto the surface using x(f" 1 ), thereby gen- 
erating mapped points on the surface; 6) sampling the 
surface using at least a subset of the mapped points to 
generate samples of the surface; and 7) analyzing the 
samples of the surface to determine characteristics of the 
surface. 

37 Claims, 14 Drawing Sheets 
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SYSTEM AND METHOD FOR ANALYZING A 
SURFACE BY MAPPING SAMPLE POINTS 
ONTO THE SURFACE AND SAMPLING THE 
SURFACE AT THE MAPPED POINTS 

FIELD OF THE INVENTION 

The present invention relates to the field of measurement 
and data acquisition, and more particularly to the analysis of 
a surface of an object by generation or mapping of a point 
set on a Riemannian manifold or surface of the object and 
analyzing data values on the mapped point set to character- 
ize the object. 

DESCRIPTION OF THE RELATED ART 

Many scientific and engineering tasks involve exploring, 
i.e., sampling or scanning, a region, such as an image, 
surface, or object, to acquire data characterizing the region. 
Examples of such tasks include parts inspection for auto- 
mated manufacturing systems, alignment tasks for auto- 
mated assembly systems, and detection, recognition, and 
location tasks in machine vision and motion control systems, 
among others. Another application relates to the determina- 
tion of integrals of smooth functions defined on Riemannian 
geometries embedded in R", i.e., in multi -dimensional real- 
space, via sums of function values at well-chosen sampling 
points. Yet another domain application relates to the explo- 
ration and characterization of unknown spaces. 

In a typical sampling or scanning system a computer 
system is coupled to a sensor, such as a camera, which is 
operable to acquire optical, image, or other information from 
a target object or surface. The computer system may take 
any of various forms. The system may also include hardware 
and corresponding software which is operable to move one 
or both of the sensor and the target object to perform the scan 
or measurement. In robotics and motion planning an under- 
standing of the underlying geometry of space is important. 
Various techniques have been developed to scan regions 
under various constraints or toward specific goals. In many 
cases the geometry is known in advance and a specific goal 
is desired, such as object recognition and/or characteriza- 
tion. In others the space in question may not be 
characterized, i.e., the nature of the space may be unknown. 
Exploration of such a space may involve goals such as 
determining the dimensionality and structure of a space with 
unknown underlying geometry, or finding critical points of 
a high-dimensional space. 

Typical applications of scanning spaces of known geom- 
etry include, but are not limited to: mine-countermeasure 
missions, continental shelf oceanographic mapping, con- 
tamination cleanup, floor scrubbing, crop plowing, non- 
destructive testing, and bridge inspections, among others. 
Many of these applications attempt efficient coverage of the 
space, and so involve determining a coverage plan for the 
scan region. Most coverage planners are still based on 
heuristics and the smoothness of the developed trajectories 
is rarely an issue. 

There are many coverage algorithms. In many cases the 
goal is to guide a robot or sensor to explore or to act within 
an environment. See, for example, J. Colgrave, A. Branch, 
"A Case Study of Autonomous Household Vacuum 
Cleaner", AIAA/NASA CIRFFSS, 1994. See also M. Ollis, 
A. Stentz, "First Results in Vision-Based Crop Line 
Tracking", IEEE International Conference on Robotics and 
Automation, 1996. 

One promising method in motion planning is based on 
Morse functions. These procedures look at the critical points 
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of a Morse function to denote the topological changes in a 
given space. See, for example, Howie Choset, Ercan Acar, 
Alfred A. Rizzi, Jonathan Luntz, "Exact Cellular Decom- 
positions in Terms of Critical Points of Morse Functions". 
5 See also Ercan U. Acar, Howie Choset, "Critical Point 
Sensing in Unknown Environments". However, Morse func- 
tions find their primary use with regard to the scanning of 
un-smooth surfaces, and so are not generally useful for many 
applications, 

10 Exploration tasks often relate to characterization of 
unknown spaces. One approach to efficient exploration of a 
surface or space is based upon Low Discrepancy Sequences 
(LDS), described below in some detail. However, current 
techniques generally involve generation of LDSs in simple 

j 5 Euclidean spaces, which may not be appropriate for more 
complex applications. 

Therefore, improved systems and methods are desired for 
generating a sampling point set, such as a Low Discrepancy 
Sequence, in a region or on a manifold or surface of an 

20 object. 

SUMMARY OF THE INVENTION 

The present invention comprises various embodiments of 
a system, method, and memory medium for analyzing a 

25 surface, wherein the surface may be a surface of an object, 
a region, and/or an n-dimensional volume. The surface may 
also be a Riemannian manifold or an n-dimensional space 
describing degrees of freedom of an object or device, such 
as a motion device, e.g., a robotic arm. The method may 

30 involve generating a point set or sequence on the surface, 
acquiring data or taking measurements of the surface at 
these sample points, and analyzing the measured data to 
characterize the surface. 

The method may include generating a sequence, such as 

35 a Low Discrepancy Sequence, within a region or space in 
which the surface is defined. The sequence may be used for 
characterizing the surface, such as a high-dimensional space, 
or for motion planning. Embodiments of the invention 
include a method for generating a Low Discrepancy 

40 Sequence in a Riemannian manifold, e.g., on an 
n-dimensional surface. 

Specifically, an embodiment of the method is described in 
which a Low Discrepancy Sequence is generated on an 
n-dimensional space, for example, a unit square, then 

45 mapped to the manifold or surface. It is also contemplated 
that other sequences of points may be mapped in a similar 
manner, and that regions other than the unit n-square may be 
used. In one embodiment, the manifold may be a Rieman- 
nian manifold. In other words, the manifold may be char- 

50 acterized by a distance-like metric. 

In one embodiment, data may be received describing a 
surface which is defined in a bounded n-dimensional space, 
also referred to as a bounded n-space. In one embodiment, 
the bounded n-dimensional space may comprise a unit 

55 square. In another embodiment, the bounded n-dimensional 
space may comprise a unit cube. As mentioned above, in 
various other embodiments, the bounded n-dimensional 
space may comprise other geometries, including unit hyper- 
cubes, or unit n-cubes of dimensionalities greater than 3, as 

60 well as n-dimensional rectangles, spheres, or any other 
geometrical region. In one embodiment, the surface may 
further be embedded in an m-dimensional space, where 
m>n. In other words, the surface may be parameterized in 
terms of a higher dimensional space according to an embed- 

65 ding function x( ). 

The received data may include a Riemannian metric 
characterizing the surface or manifold. In one embodiment, 
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the surface may comprise an image, and the received data cation relates to the location of an acoustic source via 

may further comprise an intensity function describing the discrete sampling of sound strength and direction, such as by 

image, e.g., describing pixel information of the image. a directional microphone. 

A diffeomorphism f of the bounded n-space may be One of the most straightforward applications for Low 

determined. As is well known in the art, a diffeomorphism 5 Discrepancy Sequences are calculations of integrals of 

may be considered a mapping between two spaces which smooth functions defined on Riemannian geometries 

preserves differentiability of a mapped function. embedded in R m . The replacement of integrals by sums of 

An inverse transform F 1 of the determined diffeomor- function values of well-chosen sampling points offers the 

phism may be computed. Then, a plurality of points in the same advantages as those based on low-discrepancy sets 

bounded n-space may be selected or determined. In one 10 defined in unit cubes, 

embodiment, the plurality of points may comprise a Low Another application field for Low Discrepancy Sequences 

Discrepancy Sequence. It should be noted that although is image processing, in which an intensity function of an 

selection of a Low Discrepancy Sequence of points is an image may be reinterpreted as a surface. In other words, the 

exemplary use of the method, any other point set may be surface may comprise an image. Sampling of the image may 

used as desired. 15 comprise using or measuring pixel values at pixel locations 

In response to the selection of the plurality of points, the corresponding to the mapped plurality of points. It should 

plurality of points may be mapped onto the surface via x(f wl ) also be noted that in this embodiment, sampling may not 

to generate a mapped plurality of points, e.g., a mapped Low necessarily involve measurement with a sensor, but rather 

Discrepancy Sequence. Note that the function x(f 2 ) may be may refer to calculation of an intensity value (or other 

a functional concatenation of the embedding function x( ) 20 pixel-related quantity) based on a function, e.g., an intensity 

and the inverse transform of the diffeomorphism, T 1 . function. The computed samples may then be analyzed to 

The surface may be sampled or measured using at least a determine characteristics of the image. In one embodiment, 

subset of the mapped plurality of points. In other words, tne image may comprise information related to an object, 

measurements may be made at one or more of the mapped 2S where the analyzing may be performed to determine one or 

plurality of points, thereby generating one or more samples morc characteristics of the object. 

of the surface. A Riemannian metric for the surface based on the inten- 
Finally, the generated samples of the surface may be sity function may be determined, as well as a diffeomor- 
analyzed to determine one or more characteristics of the phism representing information contained in the intensity 
surface. For example, in one embodiment, the surface may 30 function. In real applications, the integrals may be replaced 
comprise a surface of an object, and the analyzing may with sums of fiinction values. More importantly, the inverse 
comprise analyzing the samples of the surface to determine diffeomorphism in conjunction with appropriate sampling 
one or more characteristics of the object. strategies in the unit square may allow an induced sampling 
It should be noted that in some embodiments, the mapped of a & VQTi ima 5 e content. Thus, an image may be considered 
plurality of points, e.g., the mapped Low Discrepancy 35 m abstract Riemannian surface or manifold. 
Sequence, may simply be output for later use by the system It should be noted that although typical images are 
or by other systems. In other words, the sampling and/or two-dimensional images, higher dimensional images are 
analyzing may be performed by other systems, or at a later also contemplated. For example, complex multi- 
time. In one embodiment, outputting the mapped points dimensional data sets may be considered images, and the 
(Low Discrepancy Sequence) may comprise storing the 40 methods described herein may be used to explore and 
Sequence for later use. In another embodiment, outputting characterize the data sets, or objects and features contained 
the mapped Low Discrepancy Sequence may comprise therein. Examples of 3 -dimensional images include 3 d 
displaying the sequence on a display device. medical images, and 3 d seismic images, among others. It is 
Thus, by using the above-described method, a plurality of contemplated that data sets represented by even higher 
points, e.g., a Low Discrepancy Sequence, generated on a 45 dimensional images (4 d and above) may be characterized 
unit square (or other suitable geometry) may be mapped to DV various embodiments of the present invention. In other 
a surface, such as an abstract surface or manifold. It should words, any data set may be represented by an image of 
be noted thai any point set generated on the unit square (or suitable dimensionality. 

other suitable geometry) may be mapped in this way. This The intensity values at selected points may describe the 
mapped plurality of points may operate to ideally represent 50 image content in an efficient manner and may be used to 
or characterize the surface, wherein measurements taken at compress an image or to characterize its content. To a certain 
these mapped points indicate characteristics of the surface. extent, the resulting description can be regarded as a gen- 
In one embodiment, the above embodiment may be used eralized edge detector that is highly optimized for a specific 
to generate a scan path for a scanning application, such as a image. Potential applications can be found in pattern match- 
six degree of freedom alignment problem, for example. 55 ing where a given template should be characterized by a 
Other applications may include determination of scanning small amount of pixels, to speed up successive search 
sequences for one or more robotic arms, where the scan operations. The use of low-discrepancy sets offers the 
space may describe degrees of freedom for the robotic arms, advantages mentioned before. 

among others. For example, the scan space for the robotic In one embodiment, data may be received describing an 

arm may be a 6-dimensional space representing 3 positional 60 image which is defined in a bounded n-dimensional space, 

degrees of freedom, and 3 directional or orientation degrees or n-space, and which may be embedded in an 

of freedom. Alternately, the degrees of freedom may relate m-dimensional space, or m-space via an embedding function 

to the degrees of freedom of joints or other articulated x( ), where m>n. In one embodiment, the bounded n-space 

components of the arm. In one embodiment, the robotic arm may comprise a unit square (i.e., n«2). In other embodi- 

may include a directional sensor, and the sampling may 65 ments the bounded n-space may comprise a unit cube (n-3), 

include measuring a vector field at one or more of the a unit hyper-cube (n*4), or a unit n-cube where n>4. It is 

mapped plurality of points. One example of such an appli- noted that the bounded n-space may not be limited to unit 
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n -cubes, but may also comprise any other suitable geometry, include a motion control element, or actuator. The computer 

including a unit n-sphere, rectangles of various system may comprise a CPU and a memory medium, or 

dimensionalities, or any other geometry, as desired. It is also programmable logic, which is operable to store a scanning 

noted that the n-space needn't be of unit size, but may be of or analysis program that implements one of the above 

any size appropriate for the application or image under 5 methods. An embodiment of each of the above invention(s) 

consideration. may also be a software program or programs stored on a 

In one embodiment, the received data may include an memory medium, 
intensity function, as mentioned above. For example, for a 

gray-scale image, the intensity function may describe the BRIEF DESCRIPTION OF THE DRAWINGS 

gray-scale pixel value at each point in the image. Similarly 10 A better understanding of the present invention can be 

for a color image, the intensity function may provide pixel obtaincd when the following detailcd descri ption of the 

data for each point m the image. preferred embodiment is considered in conjunction with the 

A difFeomorphism (f,g) of the n-space may then be following drawings, in which: 

determined. As is well known in the art, a diffeomorphism FIGS 1B ilhlstrate rcprcscntative instrumentation 

may be considered a mapping between two spaces which and aulomation systems iocludi various measurement 

preserves differentiability of a mapped function. options- 

An inverse transform (r VJ) of the determined diffeo- FIG. 2 is a block diagram of the computer system of 

morphism may be computed. Then, a plurality of points in FIGS 1A and 1B> suitable for implementing various 

the bounded n*pace may be selected or determined. In one JQ embodiments of ^ inve ntion. 
embodiment, the plurality of points may comprise a Low 

Discrepancy Sequence. It should be noted that although **Q. 3A Ulustrates a scanning system, according to one 

selection of a Low Discrepancy Sequence of points is an embodiment of the present invention; 

exemplary use of the method, any other point set may be FIG - 3B illustrates a machine vision application of a 

used as desired. scanning system, according to one embodiment; 

After the selection of the plurality of points, the plurality FIG * 3C illuslr ates a robotics application of a scanning 

of points may be mapped onto the surface via x(r\g _1 ) to system, according to one embodiment of the present inven- 

generate a mapped plurality of points, e.g., a mapped Low *ion; 

Discrepancy Sequence. Note that the function xff^g" 1 ) PIG. 3D illustrates a phased array control application of 

may be a functional concatenation of the embedding func- 30 a scanning system, according to one embodiment of the 

tion x( )and the inverse transform of the diffeomorphism, present invention; 

(fSg -1 )- FIG. 3E illustrates an optical fiber alignment system, 

Finally, the mapped plurality of points may be analyzed to according to one embodiment of the present invention; 

determine one or more characteristics of the image. In one FIG. 4 A illustrates a Halton Sequence and a random 

embodiment, the mapped plurality of points may be ana- 35 distribution of points; 

lyzed to detect edges in the image. In other embodiments, FIG. 4B illustrates'a Low Discrepancy Sequence and an 

analyzing the mapped plurality of points may include per- Induced ^ Discrepancy Sequence, according to one 

forming pattern recognition, where features or objects com- embodiment; 

prised or represented in the image may be detected and/or A ~ ^ t „ . „ c c L . T _. 

characterized 40 illustrates Fourier Transforms of the Low Dis- 
crepancy Sequence and the Induced Low Discrepancy 

Yet another application of the method relates to a fiber Sequence of FIG. 4B, according to one embodiment; 

optic array alignment problem. Given an array of fibers CTri - A .« . . „ ,. A 4 , 

, u a * *u cu * 1. v j rIG. 5A illustrates a well-distributed point set on the 

where a second array of theses fibers must be aligned £ r . , \. 

+ i t* * ii i i -i * t i ■ . -f r surtace ot a torus, according to one embodiment; 

accurately. Typically, lasers are applied and the intensity of . 

these beams may be measured to align the arrays. 45 FIG * 5B illustrates a well-distributed point set in a unit 

0 iL c . ball, according to one embodiment; 

Said another way, the surface may comprise a multi- & ' 

dimensional space describing degrees of freedom for a FIG " 6 flowcharts a method for analyzing a surface using 

motion stage, or for first and/or second motion stages, where a ma PP ed sequence of sample points, according to one 

the first and second motion stages couple respectively to first embodiment; 

and second optical fibers (or groups of optical fibers). A laser 50 FIG - 7 illustrates edge detection in an image, according to 

beam may be directed through the first optical fiber. Samples one embodiment of the present invention; 

may be taken at a plurality of points in the search space, FIG, 8 flowcharts a method for analyzing an image using 

where the sampling includes measuring the laser beam a mapped sequence of points, according to one embodiment; 

intensity through the second optical fiber at one or more of $s and 

the mapped plurality of points. One or both of the first and 55 FIG. 9 illustrates a simplified model of an alignment 

second motion stages may be moved until the measured problem in fiber tics? according to one embo diment. 

laser beam intensity indicates that the optical fibers are . 

substantially aligned. While the invention is susceptible to various modifica- 

f . . . . , , . . uons and alternative forms, specific embodiments thereof 

It should be noted that the above examples are meant to 60 arc shown by way of e k in mc dfawi and afC ncrem 

be illustrative, and are not intended to limit the domain of described in detail It should be understood, however, that 

applicability of the method. lhc dtdLwings and detailcd 

description thereto are not 

A system may implement any of the above methods for intended to limit the invention to the particular form 

generating a Low Discrepancy Sequence in a Riemannian disclosed, but on the contrary, the intention is to cover all 

Manifold or surface and for scanning or analyzing the 65 modifications, equivalents and alternatives falling within the 

surface. The system may comprise a computer system spirit and scope of the present invention as defined by the 

coupled to a sensor. An embodiment of the system may also appended claims. 
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DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 
Incorporate by Reference 

U.S. patent application Ser. No. 09/876,982 titled "Sys- 
tem and Method for Generating a Low Discrepancy Curve 
on an Abstract Surface" which was filed on Jun. 8, 2001, and 
whose inventors are Lothar Wenzel, Ram Rajagopal, and 
Dinesh Nair, is hereby incorporated by reference in its 
entirety as though fully and completely set forth herein. 
FIGS. 1A and IB — Instrumentation and Industrial Automa- 
tion Systems 

FIGS. 1 A and IB illustrate exemplary systems which may 
store or use programs for analyzing a surface, such as a 
surface of an object, as described herein. As shown, possible 
applications of the present invention include image 
processing/machine vision, instrumentation, industrial 
automation, process control, or other purposes. These pro- 
grams may of course be stored in or used by other types of 
systems as desired. 

FIG. 1 A illustrates an instrumentation control system 100. 
The system 100 comprises a host computer 102 which 
connects to one or more instruments. The host computer 102 
may comprise a CPU, a display screen, memory, and one or 
more input devices such as a mouse or keyboard as shown. 
The computer 102 connects through the one or more instru- 
ments to analyze, measure, or control a unit under test 
(UUT) or process 150. The host computer 102 may store 
computer programs for analyzing a surface of an object as 
described herein 

The one or more instruments may include a GPIB instru- 
ment 112 and associated GPIB interface card 362, a data 
acquisition board 184 and associated signal conditioning 
circuitry 124, a VXI instrument 116, a PXI instrument 118, 
a video device 132 and associated image acquisition card 
134, a motion control device 136 and associated motion 
control interface card 138, and/or one or more computer 
based instrument cards 142, among other types of devices. 

The GPIB instrument 112 may be coupled to the computer 
102 via the GPIB interface card 362 provided by the 
computer 102. In a similar manner, the video device 132 
may be coupled to the computer 102 via the image acqui- 
sition card 134, and the motion control device 136 may be 
coupled to the computer 102 through the motion control 
interface card 138. The data acquisition board 184 may be 
coupled to the computer 102, and may interface through 
signal conditioning circuitry 124 to the UUT. The signal 
conditioning circuitry 124 may comprise an SCXI (Signal 
Conditioning extensions for Instrumentation) chassis com- 
prising one or more SCXI modules 126. 

The GPIB card 362, the image acquisition card 134, the 
motion control interface card 138, and the DAQ card 184 
may be plugged in to an I/O slot in the computer 102, such 
as a PCI bus slot, a PC Card slot, or an ISA, EISA or 
MicroChannel bus slot provided by the computer 102. 
However, these cards 362, 134, 138 and 184 are shown 
external to computer 102 for illustrative purposes. 

The VXI chassis or instrument 116 may be coupled to the 
computer 102 via a VXI bus, MXI bus, or other serial or 
parallel bus provided by the computer 102, The computer 
102 may include VXI interface logic, such as a VXI, MXI 
or GPIB interface card (not shown), which interfaces to the 
VXI chassis 116. The PXI chassis or instrument may be 
coupled to the computer 102 through the computer's PCI 
bus. 

A serial instrument (not shown) may also be coupled to 
the computer 102 through a serial port, such as an RS-232 
port, USB (Universal Serial bus) or IEEE 1394 or 1394.2 
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bus, provided by the computer 102. In typical instrumenta- 
tion control systems an instrument will not be present of 
each interface type, and in fact many systems may only have 
one or more instruments of a single interface type, such as 

5 only GPIB instruments. 

The instruments may be coupled through one or more 
sensors or transducers to the unit under test (UUT) or 
process 150. The system 100 may be used in a data acqui- 
sition and control application, in a test and measurement 

10 application, a process control application, a man-machine 
interface application, a hardware-in-the-loop applications, 
or a simulation application, among others. 

FIG. IB illustrates an exemplary industrial automation 
system 160. The industrial automation system 160 is similar 

15 to the instrumentation or test and measurement system 100 
shown in FIG. 1A. Elements which are similar or identical 
to elements in FIG. lAhave the same reference numerals for 
convenience. The system 160 comprises a computer 102 
which couples to one or more devices or instruments. The 

20 computer 102 may comprise a CPU, a display screen, 
memory, and one or more input devices such as a mouse or 
keyboard as shown. The computer 102 may connect through 
the one or more devices to a process or device 150 to 
perform an automation function, such as MMI (Man 

25 Machine Interface), SCAD A (Supervisory Control and Data 
Acquisition), portable or distributed data acquisition, pro- 
cess control, advanced analysis, or other control. In FIG. IB, 
the computer 102 may store computer programs for analyz- 
ing a surface as described herein. 

30 The one or more devices may include a data acquisition 
board 184 and associated signal conditioning circuitry 124, 
a PXI instrument 118, a video device 132 and associated 
image acquisition card 134, a motion control device 136 and 
associated motion control interface card 138, a fieldbus 

35 device 170 and associated fieldbus interface card 172, a PLC 
(Programmable Logic Controller) 176, a serial instrument 
182 and associated serial interface card 184, or a distributed 
data acquisition system, such as the Fieldpoint system 
available from National Instruments, among other types of 

40 devices. 

The DAQ card 184, the PXI chassis 118, the video device 
132, and the image acquisition card 136 may be connected 
to the computer 102 as described above. The serial instru- 
ment 182 may be coupled to the computer 102 through a 

45 serial interface card 184, or through a serial port, such as an 
RS-232 port, provided by the computer 102. The PLC 176 
may couple to the computer 102 through a serial port, 
Ethernet port, or a proprietary interface. The fieldbus inter- 
face card 172 may be comprised in the computer 102 and 

50 may interface through a fieldbus network to one or more 
fieldbus devices. Each of the DAQ card 184, the serial card 
184, the fieldbus card 172, the image acquisition card 134, 
and the motion control card 138 are typically plugged in to 
an I/O slot in the computer 102 as described above. 

55 However, these cards 184, 172, 134, and 138 are shown 
external to computer 102 for illustrative purposes. In typical 
industrial automation systems a device will not be present of 
each interface type, and in fact many systems may only have 
one or more devices of a single interface type, such as only 

60 PLCs. The devices are coupled to the device or process 150. 
FIG. 2 — Computer System Block Diagram 

FIG. 2 is a block diagram of an exemplary computer 
system 102 which may be used to implement various 
embodiments of the present invention. As FIG. 2 shows, the 

65 computer system 102 may include a memory 204 which is 
operable to store one or more software programs for imple- 
menting the methods described below. The computer system 
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102 may also include a processor or CPU 208 which may be tion of scan paths on abstract data sets see U.S. patent 

operable to execute the one or more software programs for application Ser. No. 09/876,982 titled "System and Method 

implementing the methods. The computer system 102 may for Generating a Low Discrepancy Curve on an Abstract 

include programmable logic instead of, or in addition to, the Surface" which was filed on Jun. 8, 2001, which was 

CPU 208 and memory 204. The computer system 102 may 5 incorporated by reference above. More specific applications 

store and/or execute one or more software programs which wherein the method of the present invention may be used 

perform the methods described below with reference to include robotics, and phased array control systems, as well 

FIGS. 5-10. as sampling or scanning related to image data, measurement 

The computer system 102 may also optionally include a data, acoustic data, seismic data, financial data, stock data, 

measurement or data acquisition device 206 for measuring 10 futures data, business data, scientific data, medical data, and 

or acquiring data of an object or unit under test. In other biometric data, among others. However, it is noted that the 

embodiments, e.g., where the methods described herein are present invention can be used for a plethora of applications 

used to analyze abstract data such as financial data, stock and is not limited to the above applications. Thus, the 

data, or futures data, the computer system 102 may not method of the present invention may be implemented in any 

include a measurement device 206. is of various types of devices and any of various types of 

The computer may further include an Input/Output (I/O) applications, 

interface (not shown) for communication with external FIGS. 3A-E — Example Applications of the Scanning Sys- 

systems, such as a display, keyboard, or mouse, a data tem 

acquisition device, or for communication with a network. FIGS. 3A-E illustrate various exemplary applications 

In one embodiment, the computer system 102 may 20 where various embodiments of the present invention may be 

include a display device, such as a monitor, as well as a used. However, it is noted that the invention is not limited 

chassis and one or more I/O devices, such as a keyboard to these applications, but rather may be used in any of 

and/or mouse. However, the computer system may take any various applications, 

of various forms. FIG. 3A— Image Processing System 

Referring again to FIGS. 1A, IB, and 2, the computer 25 FIG. 3A illustrates an image processing system (or scan- 
system 102 may include a memory medium on which one or ning system) according to one embodiment of the present 
more computer programs according to the present invention invention. As FIG. 3A shows, computer system 102A may 
may be stored. The term "memory medium" is intended to be coupled to a camera or sensor 132 which may be operable 
include an installation medium, e.g., a CD-ROM, floppy to acquire optical, image, or other information from an 
disks 104, or tape device; a computer system memory or 30 object 360. In one embodiment, the scanning system may 
random access memory such as DRAM, SRAM, EDO also include a motion control component 125, such as a 
RAM, Rambus RAM, etc.; or a non-volatile memory such as motion control stage, which may be operable couple to the 
a magnetic media, e.g., a hard drive, or optical storage. The computer system 102A. The motion control state may 
memory medium may comprise other types of memory as receive control signals from the computer system 102A and 
well, or combinations thereof. 35 may move the object 360 with respect to the camera/sensor 

In addition, the memory medium may be located in a first 132 for scanning purposes. In another embodiment, the 

computer in which the programs are executed, or may be motion control component may move the camera/sensor 132 

located in a second different computer which connects to the instead of the object 360 to scan the object 360. It is noted 

first computer over a network, such as the Internet. In the that as used herein, the term "scanning" may refer to 

latter instance, the second computer provides the program 40 continuous scanning or to discrete scanning, also referred to 

instructions to the first computer for execution. as sampling. 

Also, the computer system 102 may take various forms. The sensor (camera 132) may also be referred to as a 
As used herein, the term "computer system" is intended to remote scanning device. In one embodiment, the sensor 132 
encompass any of various devices which include a processor may be a camera or other image acquisition device, which 
that can execute instructions from a memory medium and/or 45 may be adapted to receive one or more portions of the 
may include a programmable logic device that can be electromagnetic (EM) spectrum, e.g., visible light, infrared, 
configured to execute a method or algorithm. Exemplary or ultraviolet (UV) light. In another embodiment, the sensor 
computer systems 102 include a personal computer system, may be an ultrasonic device for detecting sound waves. In 
mainframe computer system, workstation, network various other embodiments, the sensor may be capable of 
appliance, Internet appliance, personal digital assistant so detecting or measuring any of a variety of measurable 
(PDA), television system, embedded device, intelligent phenomena. Thus the sensor may use any of various tech- 
instrument, a smart sensor, and other systems. The computer niques to sample or scan an object, image, or surface, and 
system 102 may be configured to execute the methods produce image data, including visible light, infrared, 
described with reference to FIGS. 5-10. ultrasonic, light interferometer, and other non-contact and 

In one embodiment, the software programs which imple- 55 contact methods, or other measurement data, 

ment the methods described herein are designed for data In one embodiment, the computer system 102A may 

acquisition/generation, measurement, analysis and/or operate to "scan" previously acquired data, e.g., scan stored 

display, machine vision, simulation, or for controlling or data in a data mining application. In this instance, the system 

modeling instrumentation or industrial automation hard- may not require sensor 132 for scanning a physical object or 

ware. Example applications where the method described 60 region, as the data being analyzed has been previously 

herein may be used include inspection systems, industrial acquired and stored. 

automation or motion control systems, measurement The computer system 102 A may include a display device, 

systems, image processing or machine vision systems, such as a monitor, as well as a chassis and one or more I/O 

robotics, phased array systems and any other application devices, such as a keyboard and/or mouse. However, the 

where it is desirable to determine characteristics of an object 65 computer system 102A may take any of various forms, 

or space by measuring (or sampling or scanning) a surface In this example, one sensor 132 (a camera in this 

of the object or space. For more information on the genera- example) scanning one object 360 is shown, but in other 
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embodiments any number of sensors or cameras 132 may be late one or more sample points based upon user input and/or 

used to scan any number of objects 360. The camera 132 sensor data. In one embodiment, the sample point calcula- 

may comprise any type of camera or device operable to tions may be performed offline in a preprocessing phase. In 

acquire images of the object 360. The objects 360 may be another embodiment, part or all of the sample point calcu- 

stationary, or may be conveyed into and out of the field of 5 la tions may be performed in real time. The computer system 

view of the one or more cameras 132. 102C may be further operable to store and execute software 

The camera 132 may be operable to output a video stream programs to maneuver the one or more robotic arms 184 and 

to computer system 102, e.g., to an image acquisition respective cameras 132 to implement the calculated sam- 

(IMAQ) device comprised in the computer system 102. The pling scheme. 

computer system 102 may then analyze the images captured In one embodiment of the system shown in FIG. 3C, 

by the image acquisition board. Alternatively, the image multiple robotic arms may be used in tandem. In this case, 

acquisition board may include an on-board processor and a cooperative sampling strategy may be required which 

memory for performing a portion or all of the image coordinates the movement of each arm 184 to collectively 

analysis. As another example, the camera 132 may be a sample the object 360. 

smart camera that includes the image processing function as FIG. 3D — Phased Array Application of the Present Inven- 

well as the image acquisition function. Thus, the smart 15 tion 

camera may include a processor and memory, or program- FIG, 3D illustrates an example application of the present 

mable logic, for performing the method described herein. invention in the area of phased array control. A phased array 

More generally, a system according to an embodiment of the typically refers to a group of antennas in which the relative 

invention may be implemented in a smart sensor. phases of the respective signals feeding the antennas are 

The computer system 102A (or smart camera or other 20 varied in such a way that the effective radiation pattern of the 

smart sensor) may use or store software which analyzes the array is reinforced in a desired direction and suppressed in 

image data received from the camera 132, and controls the undesired directions. As FIG. 3D shows, computer system 

sampling or scanning operation, i.e., determines the scan 102 D may couple to a phased array 302. The phased array 

path or sample points over the object. The system may also 302 may comprise a plurality of array elements 304 which 

comprise one or more motion systems for moving the 25 mav eacn be controlled independently or in concert with the 

camera, the object, or both, in order to implement the other elements 304. The computer system 102D may 

scanning scheme. The motion system may be a motion store and execute software which is operable to control the 

control system, such as is available from National Instru- P^ 5 ^ arra y elemen Jf *> accomplish a specific task. Other 

merits Corporation examples of controlled phased arrays include telescope 

FIG. 3B — Machine Vision Application of the Present Inven- 30 ^ md ^o-mirror assemblies on fiber optic transfer 

tl °? ,. t « . ... FIG. 3E — Optical Fiber Alignment Application of the 

In a machine vision or automated inspection application present Invention 

of the present invention, shown in FIG. 3B, a system similar FIG. 3E illustrates an example machine motion applica- 

to that shown in FIG. 3A may implement the present fion where the goal is a fast procedure for precise alignment 

scanning methodology in software and/or hardware for 35 of two optical fibers In mis example application, a laser 

quality control in a manufacturing process. As FIG. 3B source 310 generates a beam 312 which is routed into a first 

shows, one or more cameras 132A-D may be coupled to fiber 3 20Aand checked or measured through a second fiber 

computer system 102B for scanmng objects from several 32OB, where the intensity of the laser beam 312 may be 

points of view. In this example, objects 360 are carried past measured or sampled at a number of points and the results 

the one or more cameras 132 by manufacturing apparatus 40 used to align the ^ fibers 320 Further details of this 

106. The system may operate to scan each object 360, and example application are provided below with reference to 

the images received from each camera 132 may be analyzed pjQ 

using image processing software executing on the computer FIGS. 4A-^C— Low Discrepancy Sequences 
system 102B. The analyses of the images may be used to Pseudo-random sequences have been used as a determin- 
detect defects or other characteristics of the object 360. For 45 ^ alternative to random sequences for use in Monte Carlo 
example, in various applications the analyses may be mcthods for differ ent problems. Recently, it was 
designed to detect one or more of: physical surface defects discovered that there is a relationship between Low Dis- 
patches, etc.); one or more components located correctly crepancy Sets and the efficient evaluation of higher- 
on the object; a correct label on the object; a correct marking dimensional integrals. Theory suggests that for midsize 
on the object; correct color information on the object, etc. 50 dimensional problems, algorithms based on Low Discrep- 
Thus, in a machine vision manufacturing application, the ancy Se ts may outperform all other existing methods by an 
results of the image analyses may be used to determine order of magnitude in terms of the number of samples 
whether an object meets desired production standards. This required to characterize the problem 
determination may be performed in any of various ways, as Given a mnction f(x> ) tne problem 0 f calculating the 
desired for a particular application. If the object does not 55 integral 
meet the desired production standards, the object may be 
rejected, as shown in the case of object 362. , 
FIG. 3C — Robotics Application of the Present Invention '(/)= 1 f(x)dx 

FIG. 3C illustrates an example application of the present 0 

invention in the field of robotics. As FIG. 2B shows, a 60 

computer system 102C may be operable to control one or in the most efficient manner is not a well posed problem. An 

more robotic arms 184, each comprising a sensor or camera approximate strategy may be based on the following proce- 

132, to scan an object 360. The computer system 102C may dure: 

be operable to store and execute software implementing a (A) Construct an infinite sequence {x 1? x^, x 3 , . . . , x, . . . } 

scanning scheme according to the present invention. More 65 of real numbers in [0, 1] that does not depend on a 

specifically, the computer system 102B may be operable to specific function f (nothing is known about f in 

store and execute one or more software programs to calcu- advance, except some general smoothness properties). 
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(B) During the n 1 * step of the algorithm calculate /(x n ) 
and the approximation to the integral in (3) as: 



(4) 



If a certain criterion is satisfied stop, else repeat step (B), 
The stopping criterion depends strongly on objectives such 
as accuracy or speed. 

This algorithm differs from standard methods such as the 
trapezoidal rule which is based on equally distributed points 
in [0, 1] in that there is no relationship between consecutive 
sets x/n)«i/n and x / (n)=i/(n+l). In other words, if the 
approximation given in equation (4) fails a design goal, a 
complete recalculation of numerous f-values is necessary. 
On the other hand, it is well known that the trapezoidal rule 
gives a 1/n 2 rate of convergence for a given continuous 
function f. 

Obviously, the quality of the trapezoidal rule is based on 
a highly homogeneous set of points. To quantify the homo- 
geneity of a finite set of points, the definition of a discrep- 
ancy of a given set was introduced: 



Z)(X) = sup|m{/?)-p(/?)| 



D{X)>B d 



(log nf d 



D(X)*B d 



(log nf 



There are also infinite sequences X with 



D(X n )$B d 



10 



(5) 



Here, R runs over all rectangles [0, r] with O^r^El, m(R) 
stands for the length r of the closed interval R, and p(R) is 
the ratio of the number of points of X in R and the number 
of all points of X. The definition given in equation (5) can 
be generalized to the case of d dimensions (d«2, 3, . . . ), 
where the term interval must be interpreted as a d dimen- 
sional rectangle. The lower the discrepancy the better or 
more homogeneous the distribution of the set. The discrep- 
ancy of an infinite sequence X = {xj, X2, ■ • , , x n , . , . ) is a 
new sequence of positive real numbers D(X„), where X„ 
stands for the first n elements of X. Other definitions for the 
discrepancy of a set exist that avoid the worst-case scenario 
according to (5). 

There exists a set of points of given length that realizes the 
lowest discrepancy. It is well known in the art that the 
following inequality (6) holds true for all finite sequences X 
of length n in the d dimensional unit cube (the Roth bound, do 
see Weisstein [1999] or Kocis and Whiten [1997]): 



25 



30 



35 



(6) 



45 



B d depends only on d. Except for the trivial case d=l, it 
is not known whether the theoretical lower bound is attain- 
able. Many schemes to build finite sequences X of length n 
do exist that deliver a slightly worse limit 



50 



(7) 



55 



(8) 



for all for all sub-sequences consisting of the first n ele- 
ments. The latter result gave rise to the definition of the Low 
Discrepancy (infinite) Sequences X. The inequality in equa- 
tion (8) must be valid for all sub-sequences of the first n 
elements, where B d is an appropriately chosen constant. 
Low Discrepancy Sequences are also known as quasi- 
random sequences. 



14 



Many of the well-studied low-discrepancy sequences in 
d-dimensional unit cubes (also referred to as unit n-cubes, 
where n is the dimensionality) can be constructed as com- 
binations of 1 -dimensional low-discrepancy sequences. The 
most popular low-discrepancy sequences are based on 
schemes introduced by Corput [1937], Halton [1960], Sobol 
[1967], and Niederreiter [1992]. 

The relationship between integrals, approximations, and 
an infinite sequence X={x 1( x^, . . . } in n dimensions is 
given by the Koksma-Hlawka inequality. 



W)-Uf)\*V{f)-D(X n ) 
/(/) = Cdxfix) 



(9) 



20 



where V(f) is the variation of f in the sense of Hardy and 
Krause. For more information regarding Low Discrepancy 
Sequences, please see the book "Random Number Genera- 
tion and Quasi-Monte Carlo Methods", by H. Niederreiter, 
CBMS-NSF Regional Conference Series in Applied Math., 
No. 63, SLAM, 1992, which gives a comprehensive intro- 
duction into the implementation of Low Discrepancy 
Sequences (Halton and Sobol'). The Halton method is 
described in detail below. All of the test results presented are 
based on Halton sequences. 

Halton sequences in 1-d start with the choice of a natural 
number greater than 1. Though not absolutely necessary, 
prime numbers p=2, 3, 5, . . . are typically chosen. If p is a 
given prime number and x„ the n** element of the Halton 
sequence, the following algorithm determines x„. 

(A) write n down in the p-ary system 

n=n^ ... no, n=n 0 +n 1 -p+ . . . n^-p* 

(B) reverse the order of the digits and add the p-ary point 
0.n 0 n! . . . n q 

(C) it is 

The n element of the Halton sequence can be calculated 
independently of all other elements. As mentioned above, in 
d dimensions one has to interpret different 1 -dimensional 
Halton sequences as coordinates of points in d dimensions. 
It is very common to start with the first d prime numbers. 

FIG. 4A shows the first 1000 elements of a Halton 
sequence (left) in the unit square (a unit n-cube of dimension 
2, i.e., n=2) for a valid choice of starting prime numbers, 
namely (2, 3). The graph on the right hand side is con- 
structed with the aid of uniformly distributed random num- 
bers. The differences are pronounced, as may be seen. 

Halton sequences in 1-D are Low Discrepancy Sets in the 
sense of equation (8). More precisely, for all n and for all 
Halton sequences X that are based on a prime number p: 



/)(X)sC p !^ with 



60 



4(p + l)logp 
P- 1 



65 



4 log p 



when p is even 
when p is odd 



(10) 



A similar result holds true for Halton sequences in d 
dimensional unit squares. In a 2-dimensional unit square for 
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the (p,q) Halton sequence with prime numbers p and q the (1.4) Let D be a well-distributed set in [0,1 ] 2 . The image of 
discrepancy is D under the transform x^^v), g _1 (u,v)) forms a well- 

distributed set lying on the aforementioned embedding in 
I (\ogn?\( p-l p+lyj-1 ?+l yi (11) R ra . 



D * X) * n * n LU log p * 2log«A2log q + 2log/tJJ S Because of (8) the latter transform preserves areas (up to 

a constant factor). In particular, (Ll)-(1.4) provides a method 

ax3 a Af* i a a t» ■ ♦ . u *• to preserve low-discrepancy properties. The algorithm (1.1)- 

FIffi.4B a nd4C-InducedPomtDistnbutions 0-4) also produces a new abstract surface: 

1 tie question arises whether an area-preserving map up to 

a constant factor from [0,1] 2 onto an abstract surface S can 1Q ^^^vj^^v^+af^^^^v^iu^G^^^^v))^ 
be found. Based on such mapping, a set of points defined in 

[0,1] 2 could be transformed onto the abstract surface, pre- where the resulting area element dA 2 is independent of u and 

serving low discrepancy in the sense of the metric defined v . Note that this expression denotes a Riemanoian metric for 

for the surface. Therefore, given a nonnegative function the abstract surface, where, as those skilled in the art know, 

*P(u,v) defined on [0,1] are there functions f(u,v), g(u,v) ^ a generalization of the concept of Euclidean distance, 

and a constant c with 0ne of the simplest examples ^ the abstract surface 

defined by 



df(u, v) dgju, v) r3/(u, v) dgju, v) (8) 

= cr(u, v) 

du dv dv du ds^dif+iPdv 1 

where (f,g) is a diffeomorphism of [0,1 ] 2 ? It is shown below 2 ° n ^™£™ ™*\<^<*^^> «fi 

that under certain conditions this choice is always possible. ^} ™ lS * mbe ^ n S * a ™ { ™ c l f , 1D P ol L ar 

Theorem l:Ut%v) be a nonnegative function on [0,1? coordinates it is %v)-u. Steps (I.l)-(1.4) deliver the 

where V 2 is continuously differentiate. Let %v) be dimorphism (f(u,v), g(u,v))=(u ,v), i.e. (f (u,v), g^u, 

positive with exception of a set L of points (u,v) of v))=(sqrt(u),v) which, of course, produces a new parameter- 

Lebesgue-measure 0. For (u, v) <£ L let 25 ization (sqrt(u)cos(2jtv), sqrt(u)sin(2jtv)) of the given circle. 

FIG. 4B compares a low-discrepancy set derived from 

1 ru f9 ) such a sequence in [0,1 ] 2 (i.e. discarding those points in 

/K V) = f 1 v)i dU * {u * V) t 0 ' 1 ? that are outside tbe circle) with the image of such a 

0 sequence under an induced transform. In both cases, a 

g{u V ) = g(v) = 1 f 1 r d udrt{u v) 30 reclan g ular re e ion in proximity to the center of an circle is 

^^dudv^u, V )Jo Jo ' depicted. The distributions look alike, however FIG. 4C 

reveals that the Fourier transforms are remarkably different; 
pronounced orientation dependence in the first case and 

Furthermore, let the functions f and g be extendable to smooth frequency behavior in the latter case. Additionally, 

continuously differentiate mappings defined on [0,1 ] 2 . 35 the first approach based on derived low-discrepancy sets in 

Then the extension functions (f,g) define a diffeomorphism [o,lf suffers from serious drawbacks in higher dimensional 

of [0,1 f/L where (8) is always valid. spaccSt For example, in high dimensional spaces, the frac- 

It is noted that the set of points L refers to the few points t ion of points outside the n-sphere approaches 1, and so the 

which are not well-behaved, and so are excluded from vast majority of the generated points are not usable, resulting 

consideration. 40 [ n substantial inefficiencies. 

Proof: The denominators m (9) are positive. The extensions A more sophisticated example is a part of a sphere given 

of f and g map [0,1] 2 onto [0,1] and are continuously as an abstract surface by 

differentiable. g(v) is strictly monotone in v (the j z m i, \ A z j z u / \ ■ 1 * j * m 

j 1 -Ft ' f\\ r* / \ t ds =cos (v)du'+d\' i where (u,v) is located in [0,1 r. 

Lebeseue -measure of L is 0). For (u, v)gL, * r»3 . jj- • / \ ) / \ ✓ v • ✓ \ / \ ■ 

e } v An R embedding is x(u,v)=(cos(u)cos(v)^in(u)cos(v),sm 

Afl , . , , 45 M). Algorithm (I.l)-{1.4) produces the diffeomorphism 

aj^oa^ _ £^£|^) . (f(u ) v),g(u > v))=(u ) sin(v)/si n (l)) ( i.e. (f-(u,v),g-(u,v)Mu, 

arcsin(sin(l)v)). Thus, a Low Discrepancy Sequence of 

d^"' v ^g'(v) - - *(m, v) points defined on the unit square may be mapped onto the 

du /oX dudv ^ v ) sphere in a manner that preserves the Low Discrepancy 

50 attributes of the sequence. 

, x , , « . v , , c c , v c/ , x r FIG- 5A — Example: Induced Point Set on a Torus 

,V e ) d imph6S V_V frDm f(U,V) " f(U V) U f ° ll0WS Another example deals with the surface of a torus. The R 3 

a ™, ° * ^'f' . — 1 1 j * i *i_ i 1 embedding and the abstract surface have the descriptions 

I neorem 1 allows the construction of well -distributed sets / . u \ 

D of points on R" that are embeddings of abstract surfaces. 55 

The procedure is based on given distributions of sets in xfcv)-(( fl+ b co S (2*v)) c 0S (2* u ), ( a+ b cos(2*v)) ^(2^5,6 

[0,1] . An R" embedding of an abstract surface (1) is a sin(2*v)) 
mapping x(u,v)«(x 1 (u,v), x 2 (u,v), . . . , x„(u,v)) with E«x w x u , 

Fo Xu x VJ and Gox^ All functions are differentiable. ds 2 -4n\a+b cosC^v^^V^Vrfv 2 

The algorithm is straightforward: 60 A1 . t . rj ^ A . . 

Algorithm I: Algorithm (UMI-4) dehvers 

0-1) Given an abstract surface S defined on [0,1 f where 

^(u^) satisfies tbe propositions of Theorem 1. In (/(«, v), g( Uj v)) = L v + ^- siu(2^v)]. 

addition, let be given an embedding x(u,v) of S in R w . ^ 1 

(1.2) According to Theorem 1 a diffeomorphism (f(u,v),g(u, 65 

v)) of [0,1 ] 2 /L is constructed. After computing the inverse transform, homogeneous 

(1.3) Compute the inverse transform (f a (u,v),g -1 (u,v)). point sets lying on the surface of a torus can be generated. 
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FIG. 5 shows such a distribution and projections onto it turns out that 
different planes. As FIG. 5 clearly shows, the induced point 
set is distributed homogeneously over the surface. 
Higher Dimensional Riemannian Geometries 

In 3- or n -dimensional Riemannian geometries the term 5 
area preserving may be replaced with volume preserving or 
content preserving. 

Volume or content preserving mappings (up to a constant 
factor) are characterized by 

■ • • .»> v <tet(g«(ui, • • ■ ,Un))=«mst. (10) 10 

It is natural to generalize Theorem 1 in the following 
manner. 

Theorem 2: Let W(u 1 ,u 2 , . . . ,u M ) be a nonnegative function 

on [0,1]" where is continuously differentiate. Let The functions f t depend only on u k and can easily be 
^(^i» u 2» ■ ■ ■ > u «) De positive with exception of a set L of 15 inverted. 

points (u 1? u 2 , . . . ,10 of Ixbesgue-measure 0. For (u 1( [t fa worthwhile t0 note ^ the alternative approach of 



/t("l,«2 «n) ! 



jH'rfUilsin*" 1 ^*)! 



^rfutlsin*- 1 ^)! 



for* = 1,2, ... ,n-l 



u^, . . . ,u„)£Llet 



/jfU!, U 2 , ... , U n ) 



/„-l(t*„-l,U n ): 



X'^Kl^Wl Un) 

Jo Jo 



rfui ... dH„_i*(«i, 



n... r'd«i ... rf«„t( W1 «„) 
Jo 



deriving a well-distributed set based on the surrounding unit 
n-cube suffers from serious limitations. The content or 
volume of such an n cube is 1, the content of an 
20 n-dimensional unit ball is 



JL 

which is on the order of 
1 



25 



30 



Furthermore, let the functions fj, be extendable to 

continuously differentiate mappings defined on [0,1]". 
Then the extension functions f=(f A ,f 2 , . . . ,f„) define a 
diffeomorphism of [0,1]7L where (10) is always valid. 
Proof: The specific structure of these functions simplifies the 
determinant of the Jacobian of f to exactly one product 

<3/i(«i « n )3/ 2 (« 2 ,... ,«„) 5/„(«„) 



du 2 



35 



The rest is similar to proof of Theorem 1. 
q.e.d. 

Algorithm I may be generalized to n dimensions as 
follows: 
Algorithm II: 

01.1) Given an abstract surface S defined on [0,1] M where 
^( u i> u 2» • • • >u„) satisfies the propositions of Theorem 2. 
Furthermore, let be given an embedding x(u l7 u 2 , . . . ,u n ) 
of S in R m , where m>n. 

(11. 2) According to Theorem 2 a diffeomorphism f(u n . . . , 

«XfiK .O, • • • ,u„» of [o,ir is 

constructed. 

(11.3) Compute the inverse transform T 1 ^, . . . ,uj. 

(11.4) Let D be a well-distributed set in [0,1]". The image of 
D under the transform x(r 1 (u lt . . . .uj) forms a well- 
distributed set lying on the aforementioned embedding in 
R"\ 

FIG. 5B — Example: Induced Point Set in a Unit Ball in 
n-Dimensional Euclidean Space 

FIG. 5B illustrates an induced point set in a unit ball 
embedded in an n-dimensional (n-3) Euclidean space. 
Given the unit ball in n-dimensional Euclidean space in 
standard parameterization 

J[u v . . . r u re )-u„(siii(2jttt 1 )3in(7u*2) * • * s^fax-i)) cos(2ju* 1 ) 
sinfjWj) . . . sinfau^) . . . ,cos(»« n _2)sin(n« J ,_ 1 ) ) cos(jtu rt _ J) 



In other words, the probability to hit the n ball is extremely 
small when n is beyond 10. Said another way, the volume in 
the unit cube which is outside the unit ball is in the corners 
of the cube. The higher the dimensionality of the unit 
n-cube, the more corners there are, thus the smaller the ratio 
of the ball volume to the n-cube volume. 

40 FIG. 6 — Flowchart of a Method for Generating a Low 
Discrepancy Sequence on a Manifold 

FIG. 6 is a flowchart one embodiment of a method for 
generating a sequence, such as a Low Discrepancy 
Sequence, on a surface or manifold, e.g., on an abstract 
45 surface. Specifically, an embodiment of the method is 
described in which a sequence or point set, e.g., a Low 
Discrepancy Sequence, is generated on an n-dimensional 
space, for example, a unit square, then mapped to the 
manifold, although it is also contemplated that other 
sequences of points may be mapped in a similar manner, and 
that regions other than the unit n-square may be used. In one 
embodiment, the manifold may be a Riemannian manifold. 
In other words, the manifold may be characterized by a 
distance-like metric. It should be noted that in various 
55 embodiments, one or more of the steps may be performed in 
a different order than shown, or may be omitted as desired. 
Furthermore, additional steps may be performed as desired. 

As FIG. 6 shows, in 602 data may be received describing 
a surface which is defined in a bounded n-dimensional 
60 space, also referred to as a bounded n-space. The data may 
comprise data values or a plurality of points on the surface. 
For example, in an image processing application, the 
received data may comprise pixel data corresponding to an 
acquired image. In a seismic signal processing application, 
65 the received data may comprise acoustic signal strengths on 
the surface or interior of a volume. In a financial application, 
the received data may comprise company stock metrics such 
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as stock price, price/earnings ratios, earnings, and volatility, words, measurements may be made at one or more of the 

among others. Similarly, the received data may include mapped plurality of points, thereby generating one or more 

futures data such as pricing, term information, and volume, samples of the surface. 

among others. In a medical application, for example, the Finally, in 614 the generated samples of the surface may 

received data may include 2- or 3-dimensional image data s be analyzed to determine one or more characteristics of the 

for x-ray charts positron emission tomography, computer- surfacc For c k ^ one embodimcnt) the 

aided I tomography, nuclear magnetic resonance , ultra-sound rise a surface of gQ obj ^ the analyzi 

or other medical images and data sets. In another medical ^ ^ { ^ ^ f ^ &urfacc 

application the received data may include Electro- / .... . . A T * . 

Cardiogram (EKG) and Electro-Encephalogram (EEG) data. ° ne ° r m °™ c ^ arac ens f ics ° f the object. In a manufacturing 

In a more abstract medical application, the received data 10 q^ahty con^ol apphcation, the object may be a circuit card, 

may include multi-dimensional diagnostic metrics. In a and , the surface of the circuit card may be sampled and 

robotics application the received data may include degrees analyzed to determine whether all of the components on the 

of freedom information regarding the operation of one or card are in the correct location and/or orientation. In a 

more robotic arms. In a computer aided identification appli- computer aided identification application the surface may be 

cation the received data may include biometric data such as 15 a human face which may be analyzed to determine an 

fingerprint data, facial structure metrics, ambulatory metrics identity of the human. In an image processing application, 

(gait information), voice data, and image data, among oth- the surface may comprise a color or gray -scale image which 

ers. In a geographic application the received data may may be analyzed for feature detection, recognition, and 

include altitude data, image data representing terrain fea- characterization. In a geographic application the surface 

tures such as temperature, moisture, reflectance, ground 20 may comprise terrain (data) which may be analyzed to 

cover, urbanization, or any other terrain feature. It is noted determine features such as altitude, temperature, moisture, 

that in each of these fields, as well as others, there are reflectance, ground cover, urbanization, or any other terrain 

applications of the present invention too numerous to men- feature. 

tion. In many applications the surface and/or object may be an 
In one embodiment, the bounded n-dimensional space 25 abstract surface or object, such as a data set describing a 
may comprise a unit square. In another embodiment, the phenomenon or relationship. For example, in a seismic 
bounded n-dimensional space may comprise a unit cube. As signal processing application, acoustic signal strengths on 
mentioned above, in various other embodiments, the the surface or interior of a volume may be analyzed to 
bounded n-dimensional space may comprise other determine the presence of geologic structures indicative of 
geometries, including unit hyper-cubes, or unit n-cubes of 30 oil or mineral deposits. In a financial application, company 
dimensionalities greater than 3, as well as n-dimensional stock metrics such as stock price, price/earnings ratios, 
rectangles, spheres, or any other geometrical region. In one earnings, and volatility, among others, may be analyzed to 
embodiment, the surface may further be embedded in an determine characteristic patterns or formations indicating 
m-dimensional space, where m>n. In other words, the sur- buy or sell conditions. Similarly, futures data such as 
face may be parameterized in terms of a higher dimensional 35 pricing, term information, and volume, among others, may 
space according to an embedding function x( ). be analyzed to determine expected values for given con- 
In one embodiment, the received data may include a tracts. 
Riemannian metric characterizing the surface or manifold. In a medical application, for example, 2- or 3-dimensional 
In another embodiment, the surface may comprise an image, image data for x-ray charts, positron emission tomography, 
and the received data may comprise an intensity function 40 computer-aided tomography, ultra-sound, or other medical 
describing the image, e.g., describing pixel information of images and data sets may be sampled and analyzed for 
the image. diagnostic purposes. In other words, specific patterns or 
In 604, a diffeomorphism f of the bounded n-space may features may be detected and characterized to diagnose a 
be determined. As is well known in the art, a diffeomorphism medical condition. For example, EKG data may be analyzed 
may be considered a mapping between two spaces which 45 to characterize the health of a patient's heart, e.g., to 
preserves differentiability of a mapped function. It is noted determine the likelihood of a heart attack. Similarly, EEG 
that the diffeomorphism f may in some embodiments be data may be analyzed to determine brain function (or 
comprised of other functions, i.e., f may be a compound or malfunction) indicating brain pathologies, such as the pres- 
concatenated function. ence of a tumor. It is also noted that non-pathological aspects 
In 606, an inverse transform f 1 of the diffeomorphism 50 of the body may be researched using various embodiments 
determined in 604 may be computed. of the present invention, such as blood flow, neurological 
Then, as indicated in 608, a plurality of points in the activity, or other physiological metrics. In other words, the 
bounded n-space may be selected or determined. In one present invention also lends itself to research, as opposed to 
embodiment, the plurality of points may comprise a Low being restricted in use for diagnostics. In a more abstract 
Discrepancy Sequence. It should be noted that although 55 medical application, multi-dimensional diagnostic metrics 
selection of a Low Discrepancy Sequence of points is an in an abstract space may be analyzed, and pattern matching 
exemplary use of the method, any other point set may be algorithms used to diagnose complex multi-system health 
used as desired. conditions or functionality in a patient. 

In 610, in response to the selection of the plurality of In a robotics application degrees of freedom information 
points in 608, the plurality of points may be mapped onto the 60 regarding the operation of one or more robotic arms may be 
surface via x(t 1 ) to generate a mapped plurality of points, analyzed to determine efficient motion paths through spaces 
e.g., a mapped Low Discrepancy Sequence. Note that the with complex constraints, such as robotically aided micro- 
function x^ 1 ) may be a functional concatenation of the surgery or tele-surgery, where a surgeon (or software 
embedding function x( ) of 602 and the inverse transform of program) at one location manipulates surgical robotic arms 
the diffeomorphism, T J of 606. 65 to perform surgery on a patient at a second location. 

In 612, the surface may be sampled or measured using at In another machine vision application, an autonomous 

least a subset of the mapped plurality of points. In other robot may map terrain information into an abstract space to 
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simplify obstacle avoidance methods. In other words, certain Direction Dependent Sensing in 3 d 

algorithms for characterizing and/or avoiding obstacles may A typical scanning problem arises when electromagnetic, 

be substantially simplified by mapping the terrain data prior acoustic, or other fields must be explored where the sensor 

to analysis. This may result in real time performance of works in a direction-selective manner and can be moved 

algorithms not feasible otherwise. 5 freely, e.g. the point of maximal strength of an acoustic field 

It should be noted that in some embodiments, the mapped in a given room must be determined and the microphone 

plurality of points, e.g., the mapped Low Discrepancy must bc directed to take a measurement. From a purely 

Sequence, may simply be output for later use by the system scanning standpoint, a 6-dimensional space must be scanned 

or by other systems. In other words, the sampling and/or efficiently. Low-discrepancy Sequences in 6 d will detect 

analyzing of 612 and 614 may be performed by other 10 mis P°"it faster on average than other methods. Moreover, 

systems, or at a later time. In one embodiment, outputting standard techniques such as grid-like scanning methods are 

the mapped points (e.g., the mapped Low Discrepancy almost unrealizable because of the unbalanced structure of 

Sequence) may comprise storing the sequence for later use. these sequences. 

In another embodiment, outputting the mapped Low Dis- Integration 

crepancy Sequence may comprise displaying the sequence is As mentioned above, one of the most straightforward 

on a display device. applications for Low Discrepancy Sequences are calcula- 

Thus, by using the above-described method, a plurality of uons of integrals of smooth functions defined on Rieman- 
points, e.g., a Low Discrepancy Sequence, generated on a nian geometries embedded in R m . The replacement of inte- 
unit square (or other suitable geometry) may be mapped to grals by sums of function values of well-chosen sampling 
a surface, such as an abstract surface or manifold. It should 20 points offers the same advantages as those based on low- 
be noted that any point set generated on the unit square (or discrepancy sets defined in unit cubes, 
other suitable geometry) may be mapped in this way. Image Processing 

In one embodiment, the surface may comprise an image, Another application field for Low Discrepancy Sequences 

and the sampling may comprise sampling the image to is image processing. Assume the intensity of an image is 

generate samples of the image. It should be noted that in this 25 described by I(u,v) where (u,v) fills the unit square, 

embodiment, sampling may not necessarily involve mea- Furthermore, let I(u,v) be sufficiently smooth in u and v. 

surement with a sensor, but rather, may refer to calculation Sucn a n image can be reinterpreted as a surface (u,v,I(u,v)) 

of an intensity value (or other pixel-related quantity) based and Theorem 1 can be applied. The Riemannian metric 

on a function, e.g., an intensity function. In another belonging to I(u,v) is 

embodiment, the sampling may involve reading pixel data 30 

from an image file. The computed samples may then be jj 2 - fi i [ <9/t "' ^ fW 2 i z 9l<M * v ^ Qi ^ u * v ^ dudv i 

analyzed to determine characteristics of the image. In one \ \ - du ) ) du dv 

embodiment, the image may comprise information related to ( fdi{u, V )\ 2> L 2 

an object, where the analyzing may be performed to deter- { \ dv ) ) 

mine one or more characteristics of the object. 35 ™ 

T t 1- ii .1* i Therefore, 

In one embodiment, the above method may be used to 

generate a scan path for a scanning application, such as a six J _ 2 v )\2 

degree of freedom alignment problem, for example, as v) = J 1 + 1 ^ J +( — ^ — j 
described below in the section titled Example Applications. 

Other applications may include determination of scanning 40 

sequences for one or more robotic arms, where the scan The resulting diffeomorphism (f(u,v),g(u,v)) represents 
space may describe degrees of freedom for the robotic arms, information contained in I(u,v). In applications, the integrals 
among others. For example, the scan space for the robotic are replaced with sums of function values. More 
arm may be a 6-dimensional space representing 3 positional importantly, the inverse diffeomorphism in conjunction with 
degrees of freedom, and 3 directional or orientation degrees 45 appropriate sampling strategies in [0,1] 2 allows an induced 
of freedom. Alternately, the degrees of freedom may relate sampling of a given image content. Thus, an image may be 
to the degrees of freedom of joints or other articulated considered an abstract Riemannian surface or manifold, 
components of the arm. In one embodiment, the robotic arm It should be noted that, although the image described 
may include a directional sensor, and the sampling may above is a two-dimensional image, higher dimensional 
include measuring a vector field at one or more of the so images are also contemplated. For example, complex multi- 
mapped plurality of points. One example of such an appli- dimensional data sets may be considered images, and the 
cation relates to the location of an acoustic source via methods described herein may be used to explore and 
discrete sampling of sound strength and direction, such as by characterize the data sets, or objects and features contained 
a directional microphone. therein. 

It should be noted that such examples are meant to be 55 Algorithm III: 
illustrative, and are not intended to limit the domain of (111.1) Compute (11) where I is a given image (e.g. gray- 
applicability of the method. levels). 

Example Applications: (III. 2) Determine the inverse diffeomorphism according to 

Exploration of an Unknown Higher Dimensional Space Theorem 1. 

In some cases the geometry of an underlying space may 60 (III. 3) Choose a Low Discrepancy Sequence in [0,1 ] 2 . 

be unknown. An example is given by a car seat which has (IH.4) Compute the image sequence of this Low Discrep- 

six degrees of freedom. In this example, some combinations ancy set under the diffeomorphism mentioned in (III. 2). 

of movement are valid and others are invalid. The goal is a The intensity values of I at these points (the mapped Low 

description of this space in terms of valid combinations (6 Discrepancy Sequence) describe the image content in an 

coordinates each). Critical points in this space may be 65 efficient manner and can be used to compress an image or to 

determined by efficient exploration of the space through the characterize its content. To a certain extent, the resulting 

use of a Low Discrepancy Sequence. description can be regarded as a generalized edge detector 
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that may be highly optimized for a specific image. Potential x(r\g _1 ) to generate a mapped plurality of points, e.g., a 

applications can be found in pattern matching where a given mapped Low Discrepancy Sequence. Note that the function 

template should be characterized by a small amount of x(r ,g _1 ) may be a functional concatenation of the embed- 

pixels, to speed up successive search operations. The use of ding function X( ) of 802 and the inverse transform of the 

low-discrepancy sets offers the advantages mentioned s diffeomorphism, (r\g _1 ) of 806. 

De fb re . Finally, in 812 the mapped plurality of points may be 

FIG. 7 demonstrates the effect of Algorithm III on a given analyzed to determine one or more characteristics of the 

image. The small rectangles represent the sampling points image. In one embodiment, the mapped plurality of points 

computed by Algorithm III. As FIG. 7 shows, the algorithm may be analyzed to detect edges in the image. In other 

identifies the edges and other significant points. Other embodiments, analyzing the mapped plurality of points may 

examples of image processing tasks suitable for various pcrforming pattern recognition, where features or 

embodiments of the present invention include pattern or objects or represented in the image may be 

objec recognition and/or detection, and object de J tec ted and/or characterized. 

characterization, among others, A t j U 0 a 1X , , . , , 

FIG. 8-Flowchart of a Method for Image Processing *f ° 0t f d ab °. Ve ' m 802 ? ™ u™ y 

FIG. 8 flowcharts one embodiment of a method for " received describing an image which is defined in a bounded 
processing and analyzing an image. It should be noted that n-dimensional space, or n-space, and which may be embed- 
as used herein, the term "image" may refer not only to a ded m an m-dimensional space, or m-space via an embed- 

2- dimensional image, but to a higher dimensional image, Action x( ), where m>n. Thus, for example, the image 
e.g., images of 3-dimensions or more, as well. Examples of ma Y , be a 2D (2-dimensional) image, and the image may be 

3 - dimensional images include 3 d medical images, and 3 d 20 considered to be embedded in a higher dimensional space 
seismic images, among others. It is also contemplated that ( 3D or higher). This higher dimensional space may be based 
data sets represented by even higher dimensional images (4 00 or derived from other characteristics or features, such as 
d and above) may be characterized by various embodiments pixel data, e.g., color and/or intensity, curvatures, gradients, 
of the present invention. In other words, substantially any abstract amplitude characteristics, Riemannian metrics, or 
data set may be represented by an image of suitable 25 otner metrics derived from the image. For example, in one 
dimensionality, where sampling may refer to extracting or embodiment, pixel data values such as intensity values, e.g., 
reading data from the data set. gray scale or color pixel data values, may be interpreted as 

As FIG. 8 shows, in 802 data may be received describing amplitude, altitude, or other features, resulting in a higher 

an image which is defined in a bounded n-dimensional dimensional space (an example of the m-dimensional space 

space, or n-space, and which may be embedded in an 30 m which the image is considered to be embedded). This 

m-dimensional space, or m-space via an embedding function higher dimensional space is suitable for various differential 

x( ), where m>n. In one embodiment, the bounded n-space geometrical methods or analysis. 

may comprise a unit square (i.e., n=2). In other embodi- The embedding strategies in image processing as 

ments the bounded n-space may comprise a unit cube (n=3), described herein enable the use of well-developed differen- 

a unit hyper-cube (n=4), or a unit n-cube where n>4. It is 35 tial geometrical methods to solve existing and/or new prob- 

noted that the bounded n-space may not be limited to unit lems. Differential geometry deals with smooth 

n-cubes, but may also comprise any other suitable geometry, (differentiable) objects in spaces of dimension 1, 2, 3, . . . 

including a unit n-sphere, rectangles of various Some of these spaces are very famUiar, e.g. n-dimensional 

dimensionalities, or any other geometry, as desired. It is also Euclidean spaces, and some of them are abstract mathemati- 

noted that the n-space is not necessarily of unit size, but may 40 cal constructions. However, even these abstract mathemati- 

be of any size appropriate for the application or image under cal spaces allow embeddings in Euclidean spaces of even 

consideration. higher dimensions. 

In one embodiment, the received data may include an In general, an embedding in R" may not be space filling, 

intensity function, as mentioned above. For example, for a This may simply mean that the object exists or lives in this 

gray-scale image, the intensity function may describe the 45 space. For example, a curve in 3D may still be considered 

gray-scale pixel value at each point in the image. Similarly, one-dimensional. However, 3 dimensions may be desirable 

for a color image, the intensity function may provide pixel °r necessary to perform the embedding. As the case of a 

data for each point in the image. It is noted that various other helic demonstrates, embedding may require at least 3 

quantities may also be represented by the intensity function. dimensions, i.e., embedding may not be possible in 2 

In 804, a diffeomorphism (f,g) of the n-space may be so dimensions (where m»2). 

determined. As is well known in the art, a diffeomorphism According to various embodiments of the invention, the 

may be considered a mapping between two spaces which field of image processing provides for countless natural and 

preserves differentiability of a mapped function. It is noted mathematically oriented embeddings. The most natural is 

that the diffeomorphism (f,g) is for illustrative purposes we interpretation of a gray-scale image as a surface in 3D. 

only, and that other representations, e.g., f, may be appro- 55 This surface may be referred to as a Monge patch. This 

priate. interpretation is based on 2 independent variables x and y 

In 806, an inverse transform (r\g -1 ) of the diffeomor- (location) and one dependent variable I(x,y)— intensity. The 

phism determined in 804 may be computed. structure: 

Then, as indicated in 808, a plurality of points in the (x, y, I(x,y)) 

bounded n-space may be selected or determined. In one 60 may form a surface in three-dimensional Euclidean space, 

embodiment, the plurality of points may comprise a Low Another example is an RGB image. The structure is 

Discrepancy Sequence. It should be noted that although (x, y, R(x,y), G(x,y), B(x,y)) 

selection of a Low Discrepancy Sequence of points is an which is a surface in 5 -dimensional Euclidean space. This 

exemplary use of the method, any other point set may be results in a 5D embedding of a color image, 

used as desired. 65 These two examples are relatively straightforward. The 

In 810, after the selection of the plurality of points in 808, following is an exemplary collection of more complex 

the plurality of points may be mapped onto the surface via embeddings. 
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(1) Given a gray-level image I(x,y), construct an artificial constraints in form of F*(x,y,z,a,p,Y)eE M^for k=l, . . . , N 
abstract surface based on the Riemannian metric: where M t are manifolds in R 6 . In case of (12) the two 

separate position and orientation spaces can (after 

ds 2 =dx*+f(x, y )dy 2 normalization) be treated separately and an induced metric 

This abstract surface may be used to generate embeddings in 5 may be constructed. Homogeneous sets in R 6 may be 

higher dimensional spaces (it is very difficult to depict these mapped onto sets with similar properties in the search space, 

abstract surfaces). Id doing so, new classes of algorithms can Topologically, the resulting search space is the product of 

be derived. two balls or spheres. Further applications of Low Discrep- 

(2) The color space of all colors can be interpreted as a Sequences in the Lie group SO(3) of all rotations in R 3 
specific Riemannian space. Based on this an embedding 10 are also contemplated. This Lie group is of fundamental 
can be constructed. importance in robotics and 3 d image processing. 

(3) Same as (1) but R, G, and B values of a color image In one embodiment, the degrees of freedom of both 
generate an unusual embedding motion stages may be represented by the search space or 

surface. In this embodiment, the space may be 

ds 2 ^R{x,y)^2G{x t y)dxdy^B{x, y )dy 2 \$ 12-dimensional, in that each motion stage has 3 positional 

(4) Instead of using intensity values, one can also start with degrees of freedom and 3 orientation degrees of freedom, 
curvature or other features. This results in a natural Networks 

embedding, where the 2D image is considered to be Tenenbaum et al. [2000], Roweis and Saul [2000] report 
embedded in a higher dimensional space (3D or higher), methods of constructing manifolds representing large data 
wherein the higher dimensional space includes one or 20 sets - ^ algorithms describe data globally where local 
more additional dimensions (additional to x and y) based properties are used. The procedures are based on nearest 
on or derived from the curvature or other features. neighbors. Included is the determination of the dimension- 
Various other types of embeddings may be used according alit y of me data under consideration. 

to various embodiments of the invention. Given a large set of discrete points in a manifold, a 

Testing and Exploration 25 well-distributed subset of points of this graph may be 

A typical scenario in nondestructive testing or during determined. Such a set should treat all parts of the point set 

search routines is based on well-chosen locations where equally well. The following algorithm can be applied. 

measurements are taken. Usually, the search space is large Algorithm IV: 

and a specific point operation is computationally or other- (I VI) Apply Theorem 2 to generate a low-discrepancy set in 

wise expensive. Efficient sampling and/or scanning strate- 30 tm s manifold. 

gies are necessary for reasonable performance. C 1 ^ 2 ) Determine nodes and edges of the graph that mini- 
It is assumed here that an efficient search strategy for an mizc me distance to this set of the manifold. To this end, 
n-dimensional cube is available. Such a strategy could be appropriate thresholds must be given, 
based on low-discrepancy sets, homogenous grids, or other ^ tne number of points of the subset under construction is 
approaches. With the aid of Theorem 2 an induced strategy 35 unknown in advance the algorithm can stop whenever a 
can be derived. certain goal is achieved, e.g. any point of the underlying 
FIG. 9— Example Application: Fiber Optic Arrays P oml set is in a neighborhood of a point constructed so far. 

FIG. 9 illustrates a simplified fiber optic array alignment Probability Distributions 

problem, given an array of fibers where a second array of Sometimes, one is interested in a well-chosen point set 

theses fibers must be aligned accurately. Typically, lasers are 40 where an underlying probability distribution is known. For 

applied and the intensity of these beams may be measured to example, find a point distribution in [0,lf where the prob- 

align the arrays. ability distribution p(u,v) is known. The connection between 

Said another way, the surface may comprise a multi- mis problem and an abstract surface is given by 

dimensional space describing degrees of freedom for first d^ 2 -d« 2 +p(u,v)W 

and second motion stages, where the first and second motion 45 

stages couple respectively to first and second optical fibers Usin S Theorem 1, appropriately chosen sets of points can 

(or groups of optical fibers). A laser beam may be directed be derived from equally distributed sequences in [0,1 ] 2 . The 

through the first optical fiber. Samples may be taken at a case of probability distributions in higher dimensional 

plurality of points in the search space, where the sampling s P aces or evcn manifolds can be treated analogously, 

includes measuring the laser beam intensity through the 50 Memory and Carrier Medium 

second optical fiber at one or more of the mapped plurality ^ e computer system 102 preferably includes a memory 

of points. One or both of the first and second motion stages medium on which software according to an embodiment of 

may be moved until the measured laser beam intensity lne P resen t invention may be stored. The memory medium 

indicates that the optical fibers are substantially aligned. ma y store one or more programs to control and execute a 

In one embodiment, if the first array is regarded as fixed 55 sc *nning operation. The memory medium may also store a 

with unknown position and orientation, the second array software program for preprocessing scan data, such as to 

may be moved to establish an alignment. A simple model of generate a Low Discrepancy Sequence, described above. In 

this scenario, shown in FIG. 9, consists of 3 position degrees one embodiment, the memory medium may also store a 

of freedom (x, y, and z) in combination with 3 orientation software program for analyzing results of a scan operation, 

degrees of freedom (pitch a, roll p, and yaw angle y). This 60 ^ c term "memory medium" is intended to include an 

results in a 6-dimensional search space (surface) where installation medium, e.g., a CD-ROM, floppy disks, or tape 

further constraints add complexity to this space. For device; a computer system memory or random access 

example, memory (RAM) such as DRAM, SRAM, EDO RAM, 

jt 2 _2<m , / \i m a \i t , v RRAM, etc.; or a non-volatile memory such as a magnetic 

*W^«nd (a-a^-fc^Y-Yc^O (12) 65 media> ^ a hard ^ 0f optical stofage ^ m ^ ofy 

where P(osilion) and O(rientation) are known parameters medium may comprise other types of memory as well, or 

and do, p 0) and y 0 are given. Amore general system may use combinations thereof. 
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In addition, the memory medium may be located in a first 
computer in which the software program is stored or 
executed, or may be located in a second different computer 
which connects to the first computer over a network, such as 
the Internet. In the latter instance, the second computer 5 
provides the program instructions to the first computer for 
execution. Also, the computer system 102 may take various 
forms, including a personal computer system, mainframe 
computer system, workstation, network appliance, Internet 
appliance, personal digital assistant (PDA), television set- 
top box, or other device. In general, the term "computer 10 
system" can be broadly defined to encompass any device 
having at least one processor which executes instructions 
from a memory medium, or any device which includes 
programmable logic that is configurable to perform a 
method or algorithm. 15 

Various embodiments further include receiving or storing 
instructions and/or data implemented in accordance with the 
foregoing description upon a carrier medium. Suitable car- 
rier media include a memory medium as described above, as 
well as signals such as electrical, electromagnetic, or digital 20 
signals, conveyed via a communication medium such as 
networks and/or a wireless link. 

However, it is noted that the present invention can be used 
for a plethora of applications and is not limited to the 
applications shown herein. In other words, the applications 25 
described herein are exemplary only, and the methods 
described herein may be used for any of various purposes 
and may be stored in and execute on any of various types of 
systems to perform any of various applications. 

Although the embodiments above have been described in 30 
considerable detail, numerous variations and modifications 
will become apparent to those skilled in the art once the 
above disclosure is fully appreciated. It is intended that the 
following claims be interpreted to embrace all such varia- 
tions and modifications. 35 

We claim: 

1 . A method for analyzing a surface, comprising: 
receiving data describing an n-dimensional surface, 

wherein the surface is defined in a bounded 
n-dimensional space, wherein the surface is embedded 40 
in an m-dimensional real space via an embedding 
function x( ), and wherein m>n; 
determining a diffeomorphism f of the n-dimensional 
space; 

computing the inverse transform f" 1 of the determined 45 

diffeomorphism f; 
selecting a plurality of points in the n-dimensional space; 
mapping the plurality of points onto the surface using 

x(f^ J ), thereby generating a mapped plurality of points 

on the surface; 
sampling the surface using at least a subset of the mapped 

plurality of points to generate samples of the surface; 

and 

analyzing the samples of the surface to determine char- S5 
acteristics of the surface. 

2. The method of claim 1, 

wherein the plurality of points comprises a Low Discrep- 
ancy Sequence; 

wherein the mapped plurality of points comprises a eo 
mapped Low Discrepancy Sequence. 

3. The method of claim 1, further comprising: 
generating output comprising the mapped plurality of 

points in response to said mapping. 

4. The method of claim 1, wherein said sampling com- 65 
prises measuring data values at the surface at a plurality of 
points in the mapped plurality of points. 



5. The method of claim 1, wherein the bounded 
n-dimensional space comprises a unit n-cube, 

6. The method of claim 5, wherein n=2, and wherein the 
unit n-cube comprises a unit square. 

7. The method of claim 5, wherein n=3, and wherein the 
unit n-cube comprises a unit cube. 

8. The method of claim 5, wherein n>3, and wherein the 
unit n-cube comprises a hyper-cube of dimensionality 
greater than 3. 

9. The method of claim 1, 

wherein the surface comprises a surface of an object; 
wherein said analyzing comprises analyzing the samples 

of the surface to determine one or more characteristics 

of the object. 

10. The method of claim 1, 

wherein the surface comprises an n-dimensional volume; 
wherein said sampling comprises sampling the volume to 

generate samples of the volume; 
wherein said analyzing comprises analyzing the samples 

of the volume to determine one or more characteristics 

of the volume. 

11. The method of claim 10, 

wherein said analyzing comprises analyzing the samples 
of the volume to determine a motion path of a sensor in 
the volume. 

12. The method of claim 1, 

wherein the surface comprises an image; 

wherein said sampling comprises sampling the image to 

generate samples of the image; 
wherein said analyzing comprises analyzing the samples 

of the image to determine characteristics of the image, 

13. The method of claim 1, 

wherein the surface comprises a Riemannian manifold; 
and 

wherein the received data comprises a Riemannian metric 
for the surface. 

14. The method of claim 1, 

wherein the surface comprises a six -dimensional space 
describing degrees of freedom of a robotic arm; and 

wherein the six dimensions correspond respectively to 
three positional degrees of freedom and three direc- 
tional degrees of freedom. 

15. The method of claim 14, 

wherein said robotic arm comprises a directional sensor; 
wherein said sampling comprises measuring a vector field 

at one or more of the mapped plurality of points with 

the directional sensor. 

16. The method of claim 1, 

wherein the surface comprises a multi -dimensional space 
describing degrees of freedom for first and second 
motion stages; 

wherein the first and second motion stages couple respec- 
tively to first and second optical fibers, wherein a laser 
beam is directed through the first optical fiber, and 
wherein said sampling comprises measuring the laser 
beam intensity through the second optical fiber at one 
or more of the mapped plurality of points, the method 
further comprising: 

moving one or both of the first and second motion 
stages until the measured laser beam intensity indi- 
cates that the optical fibers are substantially aligned. 

17. A method for analyzing a surface, comprising: 
receiving data describing a surface, wherein the surface is 

defined in a bounded n-dimensional space, wherein the 
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surface is embedded in an m -dimensional real space via 
an embedding function x( ), and wherein m>n; 
determining a diffeomorphism f of the n-dimensional 
space; 

computing the inverse transform f 1 of the determined s 
diffeomorphism f; 

selecting a Low Discrepancy Sequence in the 
n-dimensional space; 

mapping the Low Discrepancy Sequence onto the surface JQ 
using x(r 1 ), thereby generating a mapped Low Dis- 
crepancy Sequence on the surface; and 

sampling the surface using at least a subset of the mapped 
Low Discrepancy Sequence to generate samples of the 
surface; and 15 

analyzing the samples of the surface to determine char- 
acteristics of the surface. 

18. The method of claim 17, 

wherein the surface comprises a surface of an object; 
wherein said analyzing comprises analyzing the samples 20 

of the surface to determine one or more characteristics 

of the object. 

19. A method for generating a Low Discrepancy Sequence 
on an abstract surface, comprising: 

selecting an abstract surface S, wherein S is defined in a 25 
unit n-cube, wherein S is embedded in R m via an 
embedding function x( ), and wherein m>n; 

determining a diffeomorphism f of the unit n-cube; 

computing the inverse transform f 1 of the determined 
diffeomorphism f; 30 

selecting a Low Discrepancy Sequence in the unit n-cube; 

mapping the Low Discrepancy Sequence onto the embed- 
ded abstract surface S using xff^ 1 ), thereby generating 
a Low Discrepancy Sequence on the abstract surface S; ^ 
and 

generating output comprising the mapped Low Discrep- 
ancy Sequence. 

20. The method of claim 19, further comprising: 
sampling the abstract surface S using at least a subset of 40 

the mapped Low Discrepancy Sequence. 

21. The method of claim 20, wherein said sampling 
comprises making a measurement at one or more points in 
the mapped Low Discrepancy Sequence. 

22. The method of claim 19, wherein n=2, and wherein the 45 
n-cube comprises a unit square. 

23. The method of claim 19, wherein n«3, and wherein the 
n-cube comprises a unit cube. 

24. The method of claim 19, wherein n>3, and wherein the 
n-cube comprises a hyper-cube of dimensionality greater 50 
than 3. 

25. A system for analyzing a surface, comprising: 
a sensor; and 

a computer which is operable to couple to said sensor, said 
computer comprising: 55 
a CPU; 

a memory medium which is operable to store program 
instructions; and 

an input for receiving data describing an n-dimensional 
surface, wherein the surface is defined in a bounded 60 
n-dimensional space, wherein the surface is embed- 
ded in an m-dimensional real space via an embed- 
ding function x( ), and wherein m>n; 
wherein the CPU is operable to execute said program 

instructions to perform: 65 

determining a diffeomorphism f of the n-dimensional 
space; 



computing the inverse transform t 1 of the determined 

diffeomorphism f; 
selecting a plurality of points in the n-dimensional 

space; 

mapping the plurality of points onto the surface using 
x^ 1 ), thereby generating a mapped plurality of 
points on the surface; 
wherein said computer and said sensor are operable to 

perform: 

sampling the surface using at least a subset of the 
mapped plurality of points to generate samples of the 
surface; and 

wherein the CPU is further operable to execute said 
program instructions to perform: 
analyzing the samples of the surface to determine 
characteristics of the surface. 

26. The system of claim 25, 

wherein the plurality of points comprises a Low Discrep- 
ancy Sequence; 

wherein the mapped plurality of points comprises a 
mapped Low Discrepancy Sequence. 

27. The system of claim 25, wherein the CPU is further 
operable to execute said program instructions to perform: 

generating output comprising the mapped plurality of 
points. 

28. The system of claim 25, wherein the bounded 
n-dimensional space comprises a unit n-cube. 

29. The system of claim 25, 

wherein the surface comprises a surface of an object; 
wherein said analyzing comprises analyzing the samples 

of the surface to determine one or more characteristics 

of the object. 

30. A system for analyzing a surface, comprising: 
a sensor; and 

a computer which is operable to couple to said sensor, said 
computer comprising: 
a CPU; 

a memory medium which is operable to store program 
instructions; and 

an input for receiving data describing an n-dimensional 
surface, wherein the surface is defined in a bounded 
n-dimensional space, wherein the surface is embed- 
ded in an m-dimensional real space via an embed- 
ding function x( ), and wherein m>n; 
wherein the CPU is operable to execute said program 

instructions to perform: 

determining a diffeomorphism f of the n-dimensional 
space; 

computing the inverse transform f" 1 of the determined 
diffeomorphism f; 

selecting a Low Discrepancy Sequence in the 
n-dimensional space; and 

mapping the Low Discrepancy Sequence onto the sur- 
face using x(f" 1 ), thereby generating a mapped Low 
Discrepancy Sequence on the surface; 

wherein said computer and said sensor are operable to 
perform: 

sampling the surface using at least a subset of the 
mapped Low Discrepancy Sequence to generate 
samples of the surface; and 
wherein the CPU is further operable to execute said 

program instructions to perform: 

analyzing the samples of the surface to determine 
characteristics of the surface. 
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31. The system of claim 30, 

wherein the surface comprises a surface of an object; 
wherein said analyzing comprises analyzing the samples 

of the surface to determine one or more characteristics s 

of the object. 

32. A system for generating a Low Discrepancy Sequence 
on an abstract surface, comprising: 

a computer which is operable to couple to said sensor, said 
computer comprising: 10 
a CPU; 

a memory medium which is operable to store program 
instructions; and 

an input for receiving data describing an n-dimensional 
abstract surface, wherein the surface is defined in a 15 
unit n-cube, wherein the surface is embedded in R m 
via an embedding function x( ), and wherein m>n; 
wherein the CPU is operable to execute said program 

instructions to perform: 2Q 

determining a diffeomorphism f of the unit n-cube; 

computing the inverse transform f" 1 of the determined 
diffeomorphism f; 

selecting a Low Discrepancy Sequence in the unit 
n-cube; 25 

mapping the Low Discrepancy Sequence onto the 
embedded abstract surface using x(r 1 ), thereby gen- 
erating a Low Discrepancy Sequence on the abstract 
surface; and 

generating output comprising the mapped Low Dis- 30 
crepancy Sequence. 

33. A memory medium containing program instructions 
for analyzing a surface, wherein said program instructions 
are executable to perform: 

receiving data describing an n-dimensional surface, 
wherein the surface is defined in a bounded 
n-dimensional space, wherein the surface is embedded 
in an m-dimensional real space via an embedding 
function x( ), and wherein m>n; 40 

determining a diffeomorphism f of the n-dimensional 
space; 

computing the inverse transform f 1 of the determined 

diffeomorphism f; 
selecting a plurality of points in the n-dimensional space; 
mapping the plurality of points onto the surface using 

x(r J ), thereby generating a mapped plurality of points 

on the surface; 

sampling the surface using at least a sub set of the mapped 50 
plurality of points to generate samples of the surface; 
and 

analyzing the samples of the surface to determine char- 
acteristics of the surface. 



35 
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34. The method of claim 33, 

wherein the plurality of points comprises a Low Discrep- 
ancy Sequence; 

wherein the mapped plurality of points comprises a 
mapped Low Discrepancy Sequence. 

35. A memory medium containing program instructions 
for analyzing a surface, wherein said program instructions 
are executable to perform: 

receiving data describing a surface, wherein the surface is 
defined in (comprised in?) a bounded n-dimensional 
space, wherein the surface is embedded in an 
m-dimensional real space via an embedding function x( 
), and wherein m>n; 

determining a diffeomorphism f of the n-dimensional 
space; 

computing the inverse transform f" 1 of the determined 
diffeomorphism f; 

selecting a Low Discrepancy Sequence in the 
n-dimensional space; 

mapping the Low Discrepancy Sequence onto the surface 
using x(f *) thereby generating a mapped Low Dis- 
crepancy Sequence on the surface; and 

sampling the surface using at least a subset of the mapped 
Low Discrepancy Sequence to generate samples of the 
surface; and 

analyzing the samples of the surface to determine char- 
acteristics of the surface. 

36. The memory medium of claim 35, 

wherein the surface comprises a surface of an object; 
wherein said analyzing comprises analyzing the samples 

of the surface to determine one or more characteristics 

of the object. 

37. A memory medium containing program instructions 
for generating a Low Discrepancy Sequence on an abstract 
surface, wherein said program instructions are executable to 
perform: 

selecting an abstract surface S, wherein S is defined in a 
unit n-cube, wherein S is embedded in R m via an 
embedding function x( ), and wherein m>n; 

determining a diffeomorphism f of the unit n-cube; 

computing the inverse transform f 1 of the determined 
diffeomorphism f; 

selecting a Low Discrepancy Sequence in the unit n-cube; 

mapping the Low Discrepancy Sequence onto the embed- 
ded abstract surface S using x(r 1 ), thereby generating 
a Low Discrepancy Sequence on the abstract surface S; 
and 

generating output comprising the mapped Low Discrep- 
ancy Sequence. 
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